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Forecast for the remainder of the Leonid storm 
season. 
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The dust trails of comet 55P/Tempel-Tuttle lead to Leonid 
storms on Earth, threatening satellites in orbit. We present 
a new model that accounts in detail for the observed 
properties of dust tails evolved by the comet at previous 
oppositions. The prediction model shows the 1767-dust trail 
closer to Earth's orbit in 2001 than originally thought; 
increasing expected peak rates for North America observers. 
Predictions for the 2002 storms are less affected. We 
demonstrate that the observed shower profiles can be 
understood as a projection of the comet lightcurve. 
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In the orbit of short period comets, mid-infrared images taken from 
satellites show a trail of dust grains that were too large to be swept 
into the comet dust tail by radiation pressure [1]. Those grains 
represent the bulk of matter lost by a comet and are thought to be 
responsible for the recent Leonid meteor storms on Earth. Only now 
has our understanding of these storms evolved enough to use the 
observations as a probe of comet mass loss in new and unique ways. 
Unlike infrared observations, grains of different masses can be 
measured independently and at high spatial resolution. 

Dust trails are a natural consequence of the dispersion in the semi- 
major axis (Aa) of the orbits after ejection, causing some grains to 
make a wider orbit than others and return later. Recent meteor shower 
observations have demonstrated that the dust trails from each return 
of the comet are narrow and often separated, because the parent 
comet orbit is never the same. McNaught and Asher [2] and Lyytinen 
and Van Flandern [3], independently following similar work by 
Kondrat'eva and Reznikov [4], have estimated the location of those 
comet dust trails near Earth orbit by calculating the planetary 
perturbations on a single test particle that is ejected at perihelion with 
just the right difference in orbital period to end up near Earth at the 
time of a given shower. From year to year, the pattern of trails moves 
in and out of Earth's orbit, because planetary perturbations differ for 
particles at different positions along the comet orbit. From this, they 
identified the oppositions of the parent comet 55P/Tempel-Tuttle that 
were responsible for recent Leonid storms (Table I). 

Observations of the Leonid showers have now accumulated to enough 
dust frail cross sections to create an accurate comet dust trail model 
for the ejecta of 55P/Tempe!-Tuttle. This unassuming comet has an 
orbital period of about 33.25 years and is special only because it has 
passed near Earth's orbit for centuries, providing a unique historical 
record of Leonid storms. Many aspects of the physical interpretation 
of the observations are relevant to comet dust trails in general. 

We measured the 1899-dust trail cross section during the 1999 
Leonid Multi-Instrument Aircraft Campaign; a NASA and USAF 
sponsored airborne mission, while flying from Israel to the Azores [5]. 
The aircraft perspective enabled us to measure simultaneously 
intrinsically faint meteors near the zenith and intrinsically bright 
meteors near the horizon. We find that the smaller grains peak earlier 
in time and have a wider profile (Figure 1) [6]. These cross sections 
are well represented by a Lorentzian shape [7] as in Equation 1 (Figure 



2). This is the first of a series of equations (Eq. 1-9) that will describe 
the dust dispersion in a trail in three dimensions. The Zenith Hourly 
Rate (ZHR) is a commonly used measure of number influx and is 
proportional to the rate of meteors observed by a visual observer 
under clear sky conditions [8], W is the full-width at half maximum of 
the ZHR profile, while k 0 max is the time of the peak in terms of solar 
longitude X 0 , which is a measure of the Earth’s position in its orbit. 
Throughout the paper, this angle will be in Equinox J2000.0. 

At the peak of the storm, the measured influx for meteoroids of visual 
magnitude V < 6.5 (or mass 2x1 0 s g [9]) was 2.8±0.4 meteoroids 
knr 2 hr 1 [10]. This corresponds [8] to a ZHR = 4,600±700 hr' (y=i. o, 
[8]) and an impact probability of 50% for the current satellite park as 
a whole. Smaller particles must have impacted in larger numbers, but 
did not result in satellite operation anomalies [11]. There is not a 
single power law over the whole mass range, as normally assumed in 
dust trail models [1], The mass power index s = 1.64±0.05 for 
meteoroids of mass less than 2x1 0 3 g (+0 magnitude), while 
intrinsically fainter meteors have larger values, increasing to s 
=1.97±0.05 for +6 magnitude meteors of mass 5x1 0 4 g [12]. Most 
of the mass is in the larger meteoroids. At least one fireball of 4 kg 
mass was observed from Leonid MAC, while Leonids up to 5 kg are 
thought to have been responsible for impacts on the Moon during the 
crossing of the same dust trail [13]. The distribution of impact flashes 
with s = 1.6±0.1 suggests that the size distribution is not changed at 
least up to 5 kg. Integrating up to this mass, the peak influx 
corresponds to 0.070 g knr 2 hr 1 . 

Similar Lorentz-shaped profiles are found also from the mid-infrared 
brightness intensity across the dust trail of comet 22P/Kopff [14]. 
The tail of the distribution has been interpreted as a separate dust 
component from grains of different size or morphology. However, the 
meteor shower shows no apparent change of the power law size 
distribution index across the Lorentz profile. We conclude that the tail 
of the distribution appears to be dynamically related to the peak and is 
not due to a separate dust component. 

Three further dust trail cross sections were obtained in November 
2000. The 1932- and 1866- dust trail were observed using the same 
intensified video cameras from a small Cessna aircraft over Florida, 
facilitated by Bo Gustafson of the University of Florida at Gainesville 
[15]. The 1733-dust trail peaked over Europe and was observed by 
llkka Yrjola in Finland using radio forward meteor scatter to measure 


the meteor rate. In Figure 3, these results are compared to visual 
observations collected by the International Meteor Organization [16]. 

These cross sections are at appropriate distances from the calculated 
trail centers to measure the dispersion of dust in the comet orbital 
plane perpendicular to Earth's orbit. Results of Lorentz profile fits are 
summarized in Table I, which includes data from historic Leonid 
showers that originated from known trails [8], Each shower represents 
a cross section at different Aa and Ar (Figure 2), and after a different 
number of revolutions N since epoch T. The observed rate is a product 
of these three factors [2], see Eq. 2: a function f(Aa) that describes 
the initial dispersion along the orbit in terms of semi-major axis, a 
function f m - 1/N that describes the subsequent dispersion due' to 
planetary perturbations and the number of revolutions (calculated 
from the relative distance between two nearby test particles), and a 
function f(Ar) that describes the dispersion in the plane of the comet 
orbit in terms of radial heliocentric distance. ZHR 0 is the peak dust 
density in a one-revolution trail. The last two functions are derived 
iteratively here by plotting the observed width and peak intensity 
(ZHR 0 f(Ar) and f(Aa), respectively) as a function of Ar and Aa. Figure 
4 shows the result. The measured width W needs to be corrected for 
the angle e h = 18.1° at which Earth crosses the trail (Eq. 3). The 
result, W E , varies with Ar and is expected to be smallest at the trail 
center. The narrowest observed historic Leonid storms imply an 

intrinsic width of only W E ° = 0.0001 3±0. 00001 AU, or 1.9±0.2x10 4 
km. 

The variation with Ar of peak intensity (Fig. 4a) and stream width (Fig. 

4b) is skewed towards negative values of Ar for both peak intensity 
and width, with comparatively narrower width and larger peak activity 
on the sunward side of the trail. The narrowest and strongest showers 
are detected when the trail position is calculated to be just outside of 
Earth s orbit. The observed trends do not comply with a cylindrical- 
symmetric Lorentz-profile dust distribution (dashed lines in Fig. 4a and 
4b) and is not Gaussian, as assumed in earlier models [2]. 

For any given functional form, there are significant discrepancies. The 
large deviation for the 1998 encounter with the 1899-dust trail is 
understood from a perturbation by Earth in the previous return of 
1965 [3]. We now measure a trail displacement of 5r = Ar obs - Ar 031 = 
0.0031 AU from the calculated position. Other discrepancies are more 
puzzling. Especially, the 1733- and 1866-dust trail encounters in 
2000, which occurred at the same calculated Ar, but resulted in 



significantly different peak intensity and width. The agreement is not 
improved by assuming that the dust density falls off (and width 
increases) with the number of revolutions N 2 (or N) as assumed by 
Lyytinen and van Flandern [3], nor with initial Aa. The latter may sound 
surprising, because comet dust trails do show such a behavior [1,14]. 
However, unlike mid-infrared images of comet dust trails, the Leonid 
showers are always measured near perihelion. 

One important clue is that the discrepancies in peak intensity and 
width deviate in sync. When the trails are too dense, they also tend to 
be too narrow. This argues against residual effects from significant 
variations in the comet activity along the orbit, or from one return to 
the next. 

We postulate that the discrepancies are due to trail shifts 5r (and 5A 0 ) 
possibly because of the particularities of comet dust grain ejection. 
McNaught and Asher [2] assume simply ejection at perihelion in the 
direction of comet motion, while Lyytinen and van Flandern [3] assume 
no ejection but high radiation pressure forces to arrive at the same 
initial Aa. However, note that the agreement in peak time and Ar 
calculated may be fortuitous because these assumptions lead to the 
same meteoroid orbit for given Aa. 


We find a smooth variation of shifts with epoch of ejection after 
matching a symmetric profile through the variation of peak intensity 
and width with Ar. The functional form that best describes the 
dispersion of dust in the heliocentric direction is Eq. 4 (solid line in Fig. 
4a), with 5r about +0.00025 AU. The equivalent width of this 
distribution (defined as integrated profile = width*peak rate) is W f ° = 
0.00060±6 AU, or 8.9±0.9x10 4 km, a factor of three larger than the 
equivalent width of 1.57 x W E ° = 0.00020±2 AU in the perpendicular 
direction. The discrepancies from this relation are of similar magnitude 
and sign for ejections dating from the same epoch. There is a 
sinusoidal variation (Eq. 5) as a function of the year of epoch from T = 
1733 until 1932 (solid line Figure 4c). The 1733 and 1866 trails 
represent the maximum and minimum of the functional trend, thus 
explaining the relatively large differences in shower width and intensity, 
despite similar Ar 031 . 

After correcting with Eq. 5, we find that the variation of width is also 
described well by an exponential curve (Eq. 6), with about half the 
scale length. With this definition of f(Ar) (Eq. 4+5), we can plot the 
corrected peak rate as a function of the initial dispersion in semi-major 



axis to find a Lorentzian shaped f(Aa) as expressed in Eq. 7, with W a = 
0. 16±0.02 AU and ZHR 0 = 0.6±0.1 xIO 5 hr'. This function represents 
the dispersion of dust along a one-revolution dust trail of Comet 
55P/Tempel-Tuttle. Note that W a does not measure a physical 
distance, but rather a dispersion in semi-major axis. 


The offset in the peak of the curve is an expected result from radiation 
pressure effects on the grains, and should be larger for smaller 
meteoroids. Indeed, we observe a logarithmic increase (Eq. 8) of the 
mass power index (s) with Aa away from the comet position (centered 
on magnitude +3.5 meteors). 


Trail shifts can also affect the time of the peak. The peak times 
calculated [2,3] differ from the observed peak times by up to + 16- 
minutes, which is the same as the range in Eq. 5. Six of eight data 
points are fitted by Eq. 9. This completes the formalism for predicting 
future Leonid returns (Eq. 1-9). Results are in Table II. 

Implications 

Several dust trails are near Earth's orbit in November of 2001 and 
2002 (Table II). Our results argue against the large dispersion and trail 
shifts that follow from numerical models by Brown and by Gockel and 
Jehn [17]. Compared to models by McNaught and Asher [2] and 
Lyytinen and Van Flandern [3], our predictions increase the importance 
of the 1767 dust trail encounter relative to that of 1866. The 1767- 
dust trail is now expected to give the highest peak rate for Earth- 
based observers, an estimated ZHR™* = 4,200 hr’. Different solutions 
for 8r introduce an uncertainty over the range 3,000-6,900 hr 1 . The 
1866-dust trail will contribute only ZHR = 2,000-3,500 and the 1699- 
dust trail 1,300-2,500 hr 1 . However, the latter storms are slightly 
wider and both will merge into a single profile with a total fluence 1.6 
times higher. Earlier estimates [2,3] had this peak 4-10 times more 
intense. The meteors will be somewhat brighter on average than during 
the storm of 1999. Other strong showers are predicted for 2002, but 
a full Moon will illuminate this next encounter and the meteors will be 
fainter on average. No further storms are predicted until 2099. 

The observed trail shifts (-0.00025 AU) are of the same order as the 
geostationary distance (0.00028 AU). In the anti-Sun direction, for 
example, the 1767-dust trail passage in 2001 causes an equivalent 
ZHR = 11,000 hr', or about 7 particles km 2 hr' with mass > 2x1 0 s g 


at the peak. At the sunward position of geostationary orbit, the 1866- 
and 1699-dust trails peak at 6,800 and 4,500 hr', respectively. 

The Moon is positioned at a relatively large distance of 0.00258 AU. In 
2001, the most significant impacts will occur when passing the 1833- 
(ZHR = 2,800 hr') and 1800-dust trails (ZHR = 900 hr'), 2 hours 
after the Earth's passage by those trails at around 14 h and 16 h UT. 
This compares to a peak influx of about ZHR = 1,100 hr' in 1999. 
Unfortunately, the Moon will be only 3 days old. In 2002, the trails will 
remain relatively far from the Moon. 

Discussion 

The shower profiles (Figs. 1 and 3) can be understood as a projection 
of the comet's light curve. Let us assume that the dust production 
rate is proportional to the water production rate. The lightcurve of 
comet 55P/Tempel-Tuttle during the 1998 return is well described by 
[18] m r = 7.5 + 35 log r (AU) (-100 to -40 days before perihelion) 
and m r = 8.5 + 20 log r (-40 days to + 100 days after perihelion), 
with heliocentric magnitude m r = m, - 5 log d (AU). 90% of all ejection 
occurs within 60 days from perihelion. Also, the water production rate 
of comets, as observed by OH radio line observations, correlates well 
with m r , without invoking additional corrections to the OH line intensity 
or the visual magnitude: log Q H20 (r) [mol/s] = 30.74±0.02 - 

0.240±0.03 m r (magn.) [19]. 

Most of the dust ejected at heliocentric distance r will end up near 
perihelion (where Earth encounters the stream) having dispersed away 
from the comet orbit to a distance Ax perpendicular to the comet 
orbital plane: Ax = V,, 1 (r) At (r). The function At (r) is the time lapse 
from ejection until perihelion, and is readily derived from the comet 
ephemeris. By making the usual assumption that the ejection velocity 
is proportional to a power of the heliocentric distance, we can 
transform Q(r) into Q(ax) as a representation of the dust dispersion 
perpendicular to the orbital plane and, after correction for projection, 
in the path of the Earth. 

The time-independent ZHR profile width can be understood because 
each particle, to first order, will return to its point of ejection after one 
return. Thus, the width measured near perihelion reflects the 
heliocentric dependence on ejection velocities and does not necessarily 
increase with orbital period. 
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the ZHR profile, one has to invoke an increase of the election 

velocity with heliocentric distance. Within the range of comet activity 

a perfect fit is provided to the intrinsic Lorentzian shape of the dust 

density in Earths path (with W e = 0.00013 AU) for: 

(AU) V no/ = '°' 22±0 ' 05 ' °- 19±0 '°3 log M (g) + 1.27*0.05 log r 

Th e actual ejection velocity includes the comet's escape velocitv 
which is about 1.4 m/s for a comet radius R c = 1.9 km [20], Hence' 
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reported mass dependence of ejection in Eq. 10 follows from the 

explanation Tor' thf d' ^ ^ ^ ^ ^ provides a natural 
explanation for the dispersion of particles in the profile and the 

implication is that the meteoroids in the ZHR profile tails were ejected 
at relatively large heliocentric distance. 

In contrast, the measured ejection velocities are an order of maonih.H* 
smaller than predicted by the Whipple model for water vapor drag of 
cometary dust grains, modified to include adiabatic expansL 
blackbody-'imited nucleus temperature, and distributed production 
throughout the coma for ejection at perihelion, and specdically for 
particle density p - 0.7 g/cm 3 [22] and R c = 1.9 km [21J: 

Log V.,(m/s) = 1.05±0.33 - 0.167log M(g) - 0.60 log r (AU) (1,) 

fector P of^'^he P ,r d f ? r | 3 3X1 °‘ 9 me,e ° r0id iS V *= 44 m/a Within 
factor of 2). The arge tolerance reflects the various versions of Eq 11 

mat are ,n use. If the dust election velocity is proportional to the gas 

ejection velocity (Eq. 11), the result does not show the Loren, z wings 

in the observed ZHR curves. ® 

One way to reconcile the Whipple model is to consider directional 
ejection from a dust jet and only the component of the mean ejection 
ve ocity vector perpendicular to the comet orbital plane Indeed a dust 
,et was observed 1 month prior .0 the 1998 perihelion passage o 
cornet 55P/Tempel-Tuttle, with an amplitude of 25“ centered on a 
north-north-eastern direction [23], The amplitude 0, the jet motion 


suggests a hot spot at +65° N, and a rotation period of 15 33±0 02 hr 
[23]. 

The observed trail displacements (Fig. 4c) and the mass dependent 
shift in the node (Fig. 1) can be understood as an effect of such jet. 
Ejection in a northerly direction explains the negative displacement in 
node. The torque exerted by the jet will cause a precession of the spin 
axis that can qualitatively account for the observed radial displacement 
5r with a semi-period of 270±80 years (Eq. 5), and in 8X 0 with a semi- 
period of 180±20 years (Eq. 9), by changing the mean direction of 
ejection at perihelion in each return. With a nuclear axis ratio larger 
than 1.5 [ 20 ], this motion is not necessary a simple sine law, hence 
the different semi-periodicities. 

Directional ejection can account for the lack of a Lorentz wing in the 
observed f(Ar). This is because the ejection vector in the comet orbital 
plane will be mostly in the direction of comet motion at large 
heliocentric distances, while nearly perpendicular to the comet motion 
vector at perihelion. The effect is to suppress the Lorentz wings. The 
3 times higher dispersion implies that the ejection velocity at 
perihelion is V 9j = 9.1 ±1.8 m/s, still short of the Whipple speed (Eq 
11 ). 

Directional ejection has the opposite effect on the distribution of dust 
in the comet orbit f(Aa). However, ejection in the direction of motion 
can not account for the full observed dispersion with Aa. Instead, a 
dispersion in perihelion distance Aq does give the correct fall off away 
from the comet if Aq is related to a difference in semi-major axis (Aa) 
relative to that of the comet according to: 

Aa = - (1+e)' 1 G M (1-p) / (V q 2 /2 - GM (1- p) /(q±Aq)) 

+ (1+e)* 1 G M / (V q 2 /2 - GM/q) (12) 

where e is the orbit eccentricity: q = a(1-e). A good fit to the data 
(solid line in Fig. 4d) follows by plotting Q(Ax) versus <Aa>, a mean of 
the two alternative possibilities of ±Aq. For the comet velocity at 
perihelion V q = 41,600 m/s and q = 0.9766 AU, the variation in Figure 
4d is matched for Aq = 6.2±0.7 * Ax and p= 7.0±0.6 x 10 4 . The 
model predicts the decay of dust density in front of the comet, where 
no data are available. 

The parameter p in Eq. 12 is the ratio of radiation over gravitational 
forces. While most of the observed dispersion is understood in terms 


of ejection velocities, the effect of radiation pressure is to shift the 
f (^a) profile to longer Aa due to an effective decrease of the 
gravitational potential. Unlike ejection velocities, the main effect will be 
along the comet orbit. The value derived from the observed shift of 
the peak 5a is valid for a visual magnitude +3.5 Leonid meteor of initial 
mass 3x1 O' 4 g (Eq. 2). From the common equation for (3 [24], we 
conclude that the average meteoroid density is p = 0.97±0.13 g/cm 3 , 
if the radiation pressure coefficient <Q pr > = 1 and the grains are 
spherical in shape. This compares well to the estimate of p ~ 0.7 
g/cm 3 from the deceleration of a Leonid fireball [22]. 

To reconcile the observed ejection speed and its increase with 
heliocentric distance with the Whipple model, we postulate that larger 
grains fall apart in the comet coma and are the main source of the 
smaller grains. Such a scenario is not unlikely given that most of the 
mass is locked up in the larger grains. In that case, the ejection 
velocities of smaller grains reflect mostly those of the larger 
meteoroids, because gas drag is not efficient far from the nucleus 
surface. In order to explain the increasing speed with Ar, grains of 
given mass need to be derived from on average larger meteoroids 
closer to the Sun. Such an effect could occur because of increased 
thermal stresses on the grains. Indeed, the large grain mass 
distribution agrees with the value of s = 1.53±0.1 (reportedly valid 
over a wide 10‘ 12 to 10 3 kg mass range) near the nucleus of comet 
1 P/Halley and expected to reflect the dust distribution shortly after 
ejection [25]. The mass distribution for small grains is consistent 
with that expected for catastrophic fragmentation, where AN(M)*M ~ 
M +k/3 A (log M), with k=0.6 for diameters smaller than one-tenth the 
diameter of the original mass [26]. Dust fragmentation in the comet 
coma is frequently implied to account for dust distributions and comet 
dust tail striae. Our meteor observations, too, show tentative evidence 
for spatial and temporal correlations that suggest breakup more than 
one return before Earth's encounter [27], 

We now have all parameters in hand to calculate the total dust mass 
loss of 55P/Tempel-Tuttle during one return. That mass is proportional 
to the equivalent dimensions of the dust trail: W r ° by 1.57 x W E ° by 
1.57 x ((a+W,) 1 5 -a 1 s ) yr., and the peak dust density derived from: 
2HR 0 = O.SiO.lxlO 5 hr’, while ZHR = 4,600 hr'corresponds 0.070 g 
km 2 hr 1 integrated up to M = 5 kg. From this, we calculate a total dust 
mass loss for each return of comet 55P/Tempel-Tuttle of 2.6 ± 0.7 
x10 10 kg. From the observed visible magnitude light curve of comet 
55P/Tempel-Tuttle we derive a total water production loss of 
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Table 1 : Dust trail parameters from past Leonid outbursts. 


Year N* Trail* f„ 


da 


dr 


5r rf 


W„ 


W„ 


2HR„, 


ZKfL 


1999 3 1099 0.38 0.130 -0.00066 +0.00020 

1998 3 1899 0.40 0.050 +0.00440 -0.0031 

1999 4 1866 0.50 0.1 18 +0.00160 +0.00003 

2000 4 1866 0.13 0.114 +0.00077 -0.00013 

2000 2 1932 0.55 0.300 -0.00120 +0.00028 

2000 8 1733 0.27 0.064 +0.00076 +0.00039 

1969 1 1932 0.95 0.934 -0.00004 +0.00037 

1966 2 1899 0.52 0.168 -0.00013 +0.00028 

1867 1 1833 1.00 0.373 -0.00014 +0.00006 

1866 4 1733 0.37 0.059 -0.00029 +0.00051 

1833 1 1799 0.95 0.174 -0.00021 +0.00021 
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*) calculations from [2, 3] 


Table 2: Forecast for the 2001+2002 encounters. 


N Year X* Time 

(J2000) (UT) 


W 

(AU) 


FWHM 

(hr) 


2HR™“ 

(hr 1 ) 


Lyytinen Asher Brown - Time 

[3] [2| [17] [2,3] 


Nov. 17, 2001 (UT) 


1 

(1965) 

235.24 

2 

(1932) 

235.37 

3 

(1899) 

235.54 

Nov. 13, 2001 

’ (UT) 

8 

(1733) 

236.12 
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(1767) 

236. 1 19 
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236.202 
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(1633) 
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(1866) 
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Nov. 17, 2002 (UT) 

1 (1965) 235.29 

Nov. 18, 2002 (UT) 

2 (1932) 

3 (1899) 235.75 

Nov. 19, 2002 (UT) 

7 (1767) 236.615 

6 (1800) 236.710 
5 (1833) 236.709 

4 (1866) 236.871 

Nov. 19, 2006 (UT) 

2 (1932) 236.620 

Nov. 18. 2007 (UT) 
(1932) 236.109 
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Fig. 1.: 1999 Leonid storm influx profiles measured by cameras 
pointed at 37° (high) and 21° (low) elevation from the aircraft window. 
They represent masses of 5x1 0 4 g and 2x1 0 3 g, respectively. Data 
from the intensified High-Definition TV camera at 90° elevation [29] 
are shown as crosses. To facilitate comparison, the dashed line copies 
the Lorentz curve fit for the high cameras to match the peak of the 
low cameras. The activity curves are scaled to match the cumulative 
influx up to the given mass that was representative of each set of 
observations. No smoothing applied. Error bars represent the 
statistical error from the number of meteors in each interval. 



Leonid shower prediction model: 

fW 12 ) 2 

ZHR = ZHR max (k 0 -X 0 max ) 2 + ( W/2 ) 2 ( 1 ) 

ZHR"~ = ZHR 0 * f m * f(Aa) * f(Ar) (2) 

W E = 2 r tan ( W/2 ) sin (e h ) (3) 

-1450 * I Ar + 8r I 

f (Ar) = 10 (4) 

5r = Ar 0 ** 1 - Ar cal = +0.00025+0.00020* sin (2n (T-1910)/270 ) AU (5) 

-600 * I Ar + 5r I 

W E (Ar) = 1.2x10' 4 10 (6) 

f (Aa) = (W a /2) 2 / ( (Aa - 0.12±0.01) 2 + (W a /2) 2 ) (7) 

s = 2.21 + 0.41 'log (Aa) (AU) (8) 


5X o = X o 0b ‘->. o cal = -0.0001 0+0.00020*s/n (2rr(T-1 91 0)/1 80 ) AU (9) 


Fig. 2.: Diagram shows the orbit of 55P/Tempel-Tuttle and the 
definition of parameters used to decribe the location and size of the 
dust trail in the prediction model (Eq. 1-9). 





Fig. 3.: Zenith Hourly Rate curves for the 2000 encounters with the 
1932-, 1 733-, and 1866-dust trails. Black dots are results from 
intensified video cameras, while crosses are radio-MS data. Open circles 
are visual observations reported in Arlt and Gijssens [16]. The solid 
line is a fit of Lorentzian profiles. The broad dashed line is the level of 
annual shower activity in non-outburst years. The narrow dashed line 
shows the predicted rate by Lyytinen and Van Flandern [3]. 




Fig. 4a: Trail cross section along a radial direction to the Sun. Open 
symbols are observed values, dark symbols show values after 
correcting for trail shifts of Fig. 4c. 4b: As Fig. 4a, for the variation of 
shower width with heliocentric distance. 4c: Trail shifts that would fit 
the observed flux to a smooth exponential behavior (dark symbols in 
fig. 4a). The open circle shows one trail equivalent width. 4d: 

Variation of trail dust density with intrinsic semi-major axis dispersion 

(variation along the comet orbit) after correction for trail shifts in Fia 
4c. y ‘ 







